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基于 伴随 方法 的 叶片 加 工 偏差 气动 灵敏 度 分 析 
罗 佳 奇 朱 亚 路 刘 多 


(北京 大 学 工学 院 航 空 航天 工程 系 ， 北 京 100871) 


摘 要 受 加 工 精度 限制 , 叶片 加 工 外 形 与 设计 外 形 存 在 一 定 偏差 , 该 类 偏差 具有 较 强 的 随机 性 , 对 叶片 气动 性 能 造成 不 
确定 性 影响 , 加 工 叶 片 外 形 偏 差 的 不 确定 性 影响 为 几何 尺度 上 的 小 扰动 问题 本 文 首先 基于 伴随 方法 建立 了 适用 于 加 工 叶 
片 外 形 偏差 气动 灵敏 度 分 析 及 气动 性 能 偏差 预 估 模 型 . 之 后 ,计算 某 二 维 涡轮 叶片 的 气动 参数 伴随 灵敏 度 , 由 此 确定 厚度 
均匀 变化 的 测试 叶片 的 气动 性 能 偏差 , 并 与 数值 模拟 结果 进行 对 比 , 验证 伴随 方法 在 叶片 气动 不 确定 性 灵敏 度 分 析 中 的 精 
确 性 和 有 效 性 . 最后, 采用 基于 灵敏 度 分 析 的 Monte Carlo 方法 研究 满足 正 态 分 布 的 批量 加 工 叶 片 的 气动 性 能 偏差 , 确定 
了 相应 的 气动 偏差 统计 模型 ， 
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Abstract 


and the manufacturing aerodynamic shape of turbomachinery blades. 


Due to the limited machine precision, there exist deviations between the nominal 
The strong randomicity of 
manufacturing variations results in uncertainties to aerodynamic performance, which is essentially a 
small perturbation problem because the manufacturing variations are micro relative to the geometric 
parameters. In the paper, the evaluation of the resultant performance uncertainty is firstly introduced 
and completed with the adjoint sensitivity. Then the performance variations of two blades with 
uniformly thickened and thinned thickness, respectively reference to the nominal blade are predicted 
by the adjoint sensitivity, which are then compared with the computational ones to demonstrate the 
sufficient accuracy and practicability of the adjoint method in uncertainty studies. Finally, Monte 
Carlo simulations are performed on a group of blade samples, the manufacturing variations of which 
subject to normal distributions. By the adjoint method, the performance variations subject to normal 


distributions are successfully obtained. 
Key words sensitivity analysis; aerodynamic performance; adjoint method; manufacturing varia- 


tion; Monte Carlo; normal distribution 


0 引 告 


在 制造 工业 中 ， 受 加 工 精度 的 影响 ,加工 叶 三 
的 外 形 与 设计 外 形 不 可 避免 地 存在 一 定 的 偏差 。 虽 
然 每 个 外 形 设计 参数 的 加 工 偏差 较 小 ， 但 是 所 有 设 
计 参 数 对 叶片 气动 性 能 的 毕 合 影响 不 容 忽 略 . 加 工 
偏差 对 叶片 气动 性 能 的 影响 研究 早 在 上 世纪 六 、 七 
十 年 代 即 得 到 开展 9。 虽然 加 工 偏差 对 叶轮 机 械 
气动 性 能 有 重要 影响 这 一 事实 早已 得 到 科技 人 员 的 
认同 , 然而 , 受制 造 业 、 实 验 成 本 等 因 系 的 制约 , 相 
关 人 研究 在 之 后 很 长 一 段 时 间 内 部 没有 得 到 实质 性 进 
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展 . 直至 21 世纪 , 随 着 数值 计算 的 发 展 和 数值 精度 
的 提高 ， 加 工 偏差 对 叶轮 机 械 气 动 性 能 不 确定 性 影 
响 再 次 成 为 国际 性 研究 热点 

本 世纪 初 ，MIT 的 研究 人 员 B-5 基于 蒙特 卡 
洛 (Monte Carlo) 分 析 方 法 建立 了 气动 偏差 的 统计 
模型 , 成 功 地 用 数值 方法 蔡 代 了 实验 方法 . 之 后 , 德 
国 FTT 的 Lange 等 站、 DLR 的 Schell 等 中 在 气动 
外 形 扰动 统计 模型 的 基础 上 采用 三 维 计算 流体 力学 
(CFD) 软件 分 析 了 气动 外 形 扰 动 对 叶轮 机 械 定 常 、 
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非 定 常 气动 性 能 的 有 影响。 在 试验 和 数值 研究 中 ， 需 
要 首先 确定 设计 参数 小 扰动 及 相应 的 气动 性 能 变化 ， 
之 后 从 大 量 数据 样本 中 统计 分 析 设 计 参 数 小 扰动 对 
气动 性 能 变化 的 影响 。 这 种 方法 需要 大 量 样本 的 试 
验 测量 或 数值 计算 , 研究 成 本 过 高 . 

研究 设计 参数 小 扰动 对 叶片 气动 性 能 的 影响 的 
基础 与 核心 是 确定 气动 性 能 对 设计 参数 变化 的 敏感 
程度 ， 即 灵敏 度 。 早 期 的 气动 研究 中 多 采用 直接 差 
分 法 (FDM) 计算 灵敏 度 。 虽然 FDM 灵敏 度数 值 
精度 高 ， 然 而 ， 对 于 复杂 气动 问题 , FDM 灵敏 度 计 
算 效 率 低 ; 此 外 ， 扰 动 步 长 的 选取 对 灵敏 度 精 度 的 
影响 也 是 该 方法 的 关键 之 一 . 目 上 世纪 七 、 八 十 年 
代 Pironneaulsl 、Jamesonil9] 成 功 地 实现 了 基于 伴随 
方法 的 灵敏 度 计算 后 ， 基 于 伴随 灵敏 度 分 析 的 内 外 
流 气动 优化 设计 研究 得 到 了 快速 发 展 0- 区 

在 计算 伴随 灵敏 度 时 , 相对 于 气动 外 形 特征 长 
度 , 网 格 扰 动 尺 度 相 对 较 小 , 在 该 扰动 范围 内 , 伴随 
灵敏 度 基 本 不 变 . 给 定 灵敏 度 时 , 由 小 扰动 引起 的 气 
动 性 能 变化 可 以 由 小 扰动 和 灵敏 度 直 接 确 定 并 与 小 
扰动 满足 线性 关系 . 确定 了 灵敏 度 , 即 可 由 设计 参数 
小 扰动 快速 、 直 接地 确定 气动 性 能 偏差 , 偏差 估计 无 
需 额 外 的 流 场 计算 ,偏差 计算 时 间 和 样本 数量 基本 
无 关 。 此 外 ， 工 业界 往往 更 为 关心 的 是 量 产 叶 片 的 
气动 性 能 ， 这 就 要 求 对 叶片 样本 的 加 工 偏差 及 气动 
不 确定 性 进行 统计 分 析 . 根据 抽样 叶片 的 加 工 俩 差 ， 
建立 合适 的 偏差 统计 模型 ， 由 灵敏 度 即 可 快速 建立 
相应 的 气动 性 能 偏差 统计 模型 . 

本 文 首 先 简要 介绍 了 伴随 方法 的 基本 原理 ， 建 
立 面向 外 形 小 扰动 的 叶片 气动 灵敏 度 分 析 及 气动 性 
能 偏差 预 估 模型 。 对 均匀 加 厚 、 减 注 叶 片 的 气动 
性 能 进行 灵敏 度 分 析 及 偏差 预 估 ， 通 过 与 数值 模 
拟 结果 的 对 比 验证 所 发 展 的 基于 伴随 方法 的 加 工 偏 
差 气 动 不 确定 性 灵敏 度 分 析 的 有 效 性 。 对 偏差 满 
足 正 态 分 布 的 批量 加 工 叶片， 采用 灵敏 度 分 析 进 
行 Monte Carlo 模拟 ， 建 立 气 动 性 能 偏差 的 统计 
模型 . 


1 伴随 方法 基本 原理 


伴随 方法 最 早 由 Pironneaulsl 、Jameson [9] 引入 
到 流体 力学 优化 设计 领域 . 由 于 其 灵敏 度 计 算 的 高 
效 性 和 高 精度 ,在 内 外 流 气 动 优化 设计 中 得 到 了 快 
速 发 展 . 伴随 方法 的 基本 原理 是 将 流动 控制 微分 方 
程 作为 约束 ,采用 伴随 算 子 将 流动 方程 的 变 分 引入 
到 目标 函数 变 分 中 ,通过 消除 流 场 变 分 对 目标 函数 
变 分 的 影响 ， 确 定 伴 随 方程 及 相应 的 灵敏 度 .。 定义 


目标 函数 7 
T= 1(w,F) (1) 


其 中 : w 为 流 场 变量 , F 为 几何 边界 . 目标 函数 的 变 
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流动 控制 方程 RE 满足 

R(w, F)=0 (3) 
其 变 分 为 

oR oR 


引入 伴随 算 子 亚 , 由 式 (2) 和 式 (4) 有 


of oR oT oR 
6T = { 药 wo} swt | 过 -wr 
(5) 
右 令 式 (5) 右 端 第 一 项 系数 为 零 , 即 可 消除 流 场 变 
分 对 目标 函数 变 分 的 影响 ， 由 此 确定 目标 函数 对 应 
的 伴随 方程 及 灵敏 度 , 分 别 如 式 (6)、 式 (7) 所 示 ， 
of OR 


Bo (6) 
0 2 TOR (7) 
6F OF oF 


由 式 (7) 可 知 : 确定 流 场 和 伴随 场 后 , 只 需 扰动 
几何 边界 实现 网 格 变 分 ， 即 可 计算 灵敏 度 ， 相 对 于 
流 场 、 伴 随 场 的 求解 ， 网 格 变 分 的 时 间 几 乎 可 以 忽 
略 不 计 . 因此 ,对 于 每 个 目标 函数 ， 只 需求 解 一 次 流 
动 控制 方程 和 一 次 伴随 方程 ， 即 可 确定 所 有 灵敏 度 
信息 , 灵敏度 计 算 时 间 和 设计 参数 的 个 数 基本 无 关 . 
在 前 期 工作 中 0" 习作 者 已 经 成 功 实现 了 伴随 灵敏 
度 的 求解 及 在 气动 优化 设计 中 的 应 用 ， 本 文 将 在 此 
基础 上 开展 研究 . 


2 加 工 偏差 灵敏 度 分 析 


加 工 叶片 外 形 偏差 的 产生 属于 随机 过 程 ， 需 要 
对 一 定数 量 的 样本 进行 特征 分 析 以 确定 加 工 系统 的 
偏差 模 态 加 , 由 此 描述 批量 加 工 叶片 的 外 形 偏差 . 由 
于 缺乏 加 工 叶片 的 外 形 数据 ， 研 究 中 仅 对 两 个 测试 
叶片 进行 气动 灵敏 度 分 析 . 在 设计 叶片 表面 均匀 分 
布 49 个 测量 点 , 并 施加 均匀 的 正 、 负 厚度 变化 生成 
测试 叶片 ,厚度 变化 为 0.2 mm, 图 1 给 出 了 某 二 维 
涡轮 叶片 的 设计 外 形 ， 以 及 均匀 加 厚 、 均 匀 减 薄 的 
测试 叶片 气动 外 形 。 为 了 减 小 叶片 尾 缘 影响 ,偏差 
测量 区 间 为 [0, 0.99cl, ec 为 叶片 驶 长 。 
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图 1 
Fig. 1 Profiles of nominal and manufacturing blades 
早期 研究 5422 发 现 : 当 外 形 扰动 量 较 小 时 , 灵 
敏 度 在 很 客 的 扰动 范围 内 基本 保持 不 变 。 加 工 叶 片 
的 外 形 偶 差 一 般 为 萤 米 级 甚至 更 低 . 因而 ， 由 加 工 
偏差 引起 的 气动 不 确定 性 分 析 属 于 小 扰动 问题 。 由 
式 (7) 可 知 ， 当 扰动 较 小 时 , 灵敏 度 近 似 为 弟 数 , 气 
动 性 能 参数 的 变化 可 以 确定 为 : 


AI= > giEi + O(e?) (8) 


其 中 : 9 为 一 阶 灵 敏 度 , < 为 外 形 偏 差 , A7 表示 气动 

为 了 验证 基于 伴随 方法 的 气动 性 能 偏差 灵敏 度 
分 析 的 精确 性 有 效 性 ， 人 研究 中 计算 了 不 同 目标 函数 
对 49 个 测量 点 的 灵敏 度 , 目标 函数 分 别 选 取 为 出 口 
流量 、 总 压 、 等 入 效率 , 定义 如 式 (9) 所 示 : 


. 1 1 

m= | dm r= 寺 Ardm 二 | 人 Amam (9) 
其 中 : dm、Ar、A7 分 别 表示 出 口 网 格 单元 的 流量 、 
压 比 和 等 精 效 率 ,， 定义 如 下 式 ; 


1 — (p/po,exit) /7 
| 三 (p/poan) -DY 


人 区 一 ce 
DPo,in 


An = (10) 
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20 
Control points 
图 2 伴随 灵敏 度 (a) 流量 ; (b) 总 压 ; (c) 等 箭 效率 


Fig. 2 Adjoint sensitivity of: (a) mass flow; (b) total pressure; (c) adiabatic efficiency 
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其 中 : p、po 分 别 为 静 压 和 总 压 ，7 为 比 热 比 , 下 标 
in、exit 分 别 表示 进 、 出 口 。 

流 场 通过 求解 Euler 方程 确定 , 进口 总 压 为 一 个 
标准 大 气压 , 总 温 为 288.15K, 出 口 等 精 马 赫 数 约 为 
0.8， 叶 片 通道 内 有 一 道 强 激 波 . 图 2 比较 了 伴随 灵 
敏 度 (AM) 和 FDM 灵敏 度 , 图 中 PS.、S.S.、LE 分 
别 表示 压力 面 、 吸 力 面 及 叶片 前 缘 . 由 图 可 知 : 对 于 
相差 一 个 量 级 的 外 形 扰动 , 灵敏 度 几乎 完全 重合 , 收 
敛 性 较 好 。 对 于 流量 目标 函数 ,伴随 灵敏 度 和 FDM 
灵敏 度 易 合 较 好 ; 对 于 压 比 和 等 粹 效率 ， 伴 随 灵敏 
度 仅 在 吸力 面 激 波 附近 与 FDM 灵敏 度 有 一 定 的 偏 
差 , 早期 研究 55 也 发 现 了 该 现象 : 受 离散 格式 、 网 
格 分 辨 率 等 因素 的 影响 ， 在 跨 超 音速 区 域 ， 基 于 连 
续 伴 随 方 法 的 灵敏 度 和 FDM 灵敏 度 会 存在 一 定 的 
偏差 . 此 外 , 在 吸力 面 尾 缘 附 近 , 灵敏 度 出 现 较 大 幅 
度 的 振荡 ， 主 要 原因 为 叶片 气动 性 能 变化 对 叶片 尾 
缘 外 形变 化 较为 敏感 . 整体 上 , 总 压 、 等 粹 效率 对 气 
动 外 形变 化 的 敏感 程度 远 高 于 流量 ， 

表 1 比较 了 三 个 不 同 叶 片 的 流量 、 总 压 和 等 炉 
效率 ， 其 中 : 设计 叶片 的 气动 参数 均 已 单位 化 ， 其 
它 叶 片 的 气动 参数 均 为 相对 于 设计 叶片 的 偏差 百 分 
比 . 在 叶片 喉 道 流 动 为 音速 ， 吃 道 面 积 直接 决定 了 
进口 流量 。 相 对 于 设计 叶片 ， 均 匀 加 厚 的 叶片 通道 
喉 道 面积 减 小 ， 导 致 流量 和 激 波 前 压力 降低 ， 流 动 
损失 增 大 ， 出 口 总 压 降 低 。 对 于 均匀 减 薄 叶片 ， 其 
通道 面积 加 大 ， 结 论 相 反 。 由 于 流量 的 伴随 灵敏 度 
精度 较 高 ， 相 对 于 设计 叶片 ， 由 灵敏 度 分 析 确 定 的 
气动 性 能 偏差 与 精确 CFD 求解 的 气动 性 能 偏差 之 
间 的 相对 误差 较 小 ， 均 不 超过 1%; 由 寺 总 压 、 效 
率 的 伴随 灵敏 度 精度 相对 较 低 ， 由 灵敏 度 分 析 确 定 
的 气动 偏差 与 精确 CFD 偏差 之 间 的 相对 误差 较为 
明显 ， 
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表 1 基于 灵敏 度 分 析 的 气动 性 能 偏差 


Table 1 Deviations of adjoint-based performance 


evaluation 
气动 参数 。 ”叶片 CFD 伴随 方法 偏差 /% 
设计 1.0000 二 一 
流量 加 厚 一 0.6336 一 0.6371 0.5489 
减 薄 0.6262 0.6239 一 0.3752 
设计 1.0000 一 四 
总 压 加 厚 一 0.1394 一 0.1489 6.8456 
减 薄 0.2068 0.1974 一 4.5385 
设计 1.0000 一 
等 炉 效 率 “加 厚 -4.8154 —4.4831 一 6.9013 
减 薄 7.1334 7.4642 4.6378 


3 Monte Carlo 统计 分 析 


Monte Carlo 方法 在 很 多 学 科 中 已 得 到 了 广泛 
应 用 。 早 期 研究 均 采 用 CFD 数值 模拟 方法 来 计 
算 抽样 叶片 的 气动 性 能 参数 并 构造 合理 的 统计 模 
型 。Garzon 等 中 采用 2000 个 样本 建立 了 气动 偏 
差 统 计 模 型 . 很 显然 ， 通 过 灵敏 度 分 析 计 算 气 动 性 
能 偏差 ,在 进行 Monte Carlo 统计 分 析 时 不 需要 额 
外 的 流 场 计算 ， 可 以 显著 提高 气动 不 确定 性 统计 分 
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析 的 效率 。 
研究 中 采用 拉丁 超 立 方 生成 满足 正 态 分 布 的 样 
本 叶片 , 样本 数 为 5000, 正 态 分 布 函 数 如 下 式 所 示 ， 


1 (cg—hg) 
e ss , és €|—E,E 
fleg)= 4 V2rog ed | (11) 


0， otherwilse 


其 中 : =0.2 mm 为 最 大 偏差 , js 和 og 分 别 为 均值 
和 均 方差 . 图 3 给 出 了 不 同 均值 的 加 工 偏差 概率 密 
度 函 数 (PDF) 分 布 , 均 方差 均 为 oo=0.375E，。 

对 抽样 产生 的 批量 叶片 , 由 式 (8) 可 以 快速 计算 
确定 气动 参数 偏差 样本 . 图 4~6 分 别 给 出 了 流量 、 总 
压 、 效 率 偏差 的 PDF 分 布 . 由 图 可 知 : 当 叶 片 加 工 
偏差 满足 正 态 分 布 时 ,气动 参 数 的 偏差 也 近似 满足 
正 态 分 布 , 如 式 (12) 所 示 ， 


1 (ef 一 pf)2 
人 or Ere 一 0 min, 0 max 
fler)= 2Aor Pe ja 


0， otherwilse 


12) 


其 中 : M、cof 分 别 为 气动 偏差 统计 模型 的 均值 和 均 
方差 ， Omin 和 Omax 分 别 为 最 小 、 最 大 偏差 ， 
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图 3 不 同 均值 的 外 形 偏差 概率 密度 函数 : (a) Hg=0.0; (b) jg= 0.58; (c) Hg = 一 0.5E 
Fig. 3 PDFs with different geometric mean values: (a) Wg=0.0; (b) ug= 0.5E; (c) Hg = —0.5E 
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图 4 流量 偏差 概率 密度 函数 ; (a) Hg= 0.0; (b) Hg= 0.5E; (c) Hg = —0.5E 
Fig. 4 PDFs of statistical mass flow deviations: (a) ug= 0.0; (b) Hg= 0.5E; (c) Hg = —0.5E 
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图 5 总 压 偏差 概率 密度 函数 : (a) Hg= 0.0; (b) jg= 0.5E; (c) Hg = 一 0.5 妃 
Fig. 5 PDFs of statistical total pressure deviations: (a) ug= 0.0; (b) ug= 0.5E; (c) Hg = —0.5E 
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图 6 等 粹 效率 偏差 概率 和 密度 函数 : (a) jg= 0.0; (b) jg= 0.5E; (c) jg = 0.5B 
Fig. 6 PDFs of statistical adiabatic efficiency deviations: (a) ug= 0.0; (b) ug= 0.5E; (c) Hg = —0.5E 


由 图 4~6 可 知 : 当 平 均 加 工 偏差 为 0.1 mm 
(ug = 土 0.5B) 时 , 与 设计 叶片 相 比 , 流量 、 总 压 、 效 
率 偏差 的 统计 平均 值 分 别 为 0.63%、0.13%、0.16%， 
此 外 ,统计 平均 的 气动 性 能 偏差 与 加 工 偏差 呈 线 性 
关系 ， 主 要 原因 是 在 本 文 的 研究 中 ， 和 气动 性 能 偏差 
只 由 一 阶 灵 敏 度 计算 确定 ， 如 式 (8) 所 示 ， 

偏 度 (Skewness) 是 统计 模型 的 另 一 个 重要 参数 ， 
其 定义 如 式 (13) 所 示 : 


_1 Ef,i — Mf . 
Skewness = 2 (于 (13) 
其 中 : n 为 样本 叶片 数 . 表 2 给 出 了 流量 、 总 压 、 效 
率 偏差 统计 模型 的 偏 度 。 由 表 2 可 知 : 最 大 偏 度 不 
超过 0.01， 进 一 步 证 实 了 气动 性 能 偏差 统计 模型 接 
近 于 正 态 分 布 。 


表 2 气动 性 能 偏差 统计 模型 的 偏 度 


Table 2 Skewness of statistical performance 


流 基 总 压 效率 
0.002377 一 0.008859 一 0.008861 


均 方 差 可 以 用 于 衡量 样本 在 抽样 空间 内 偏离 平 
均 的 程度 。 文 中 还 研究 了 三 种 不 同 均 方差 正 态 分 
布 的 外 形 偏差 统计 模型 所 对 应 的 气动 偏差 统计 模 


型 的 均 方差 变化 ， 如 表 3 所 示 。 三 种 正 态 分 布 外 
形 偏差 统计 模型 的 均值 均 为 零 , 均 方差 分 别 为 0.25 友 、 
0.375 马 、0.50。 由 表 可 知 : 随 着 外 形 偏差 统计 模型 
均 方 差 的 增 大 ， 气 动 偏 差 统计 模型 的 均 方差 也 随 着 
增 大 . 当 外 形 偏差 集中 在 均值 附近 时 ， 样 本 叶片 的 
气动 性 能 和 设计 叶片 之 间 的 差别 较 小 ， 即 气动 性 能 
也 集中 在 设计 叶片 附近 ， 从 而 使 得 均 方 差 较 小 。 上 
述 结果 也 较 好 地 反映 了 气动 偏差 统计 模型 与 外 形 偏 
差 统 计 模型 之 间 的 线性 关系 ， 


表 3 不 同 外 形 标准 差 的 气动 偏差 统计 模型 的 均 方差 
Table 3 Standard deviations of the statistical 


performance with different geometric variances 


气动 参数 o=0.25E oc=0.375E o=0.50E 
流量 /10 一 " 1.223 1.716 2.446 
总 压 /10-3 0.985 1.519 1.970 
效率 /10-3 1.160 1.857 2.320 
> 
4 结 论 


本 文 主要 介绍 了 基于 伴随 方法 的 叶片 气动 不 确 
定性 灵敏 度 分 析 方 法 并 采用 Monte Carlo 方法 建立 
了 气动 性 能 偏差 统计 模型 。 通 过 数值 模拟 及 基于 伴 
随 方法 的 灵敏 度 分 析 ， 确 定 了 均匀 加 厚 、 均 匀 减 薄 
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两 个 测试 叶片 的 气动 性 能 参数 ， 并 与 设计 叶片 进行 
了 对 比 。 最 后 开展 了 基于 灵敏 度 分 析 的 Monte Carlo 
模拟 , 分 析 了 满足 不 同 正 态 分 布 的 外 形 偏差 统计 模 
型 的 气动 性 能 参数 变化 ， 并 确定 了 相应 的 气动 偏差 
统计 模型 . 主要 结论 有 : 

1) 相对 于 总 压 、 效 率 目标 函数 , 流量 目标 郴 数 
的 伴随 灵敏 度 精 度 更 高 ， 因 而 基于 灵敏 分 析 预 测 得 
叶片 流量 偏差 精度 也 较 高 ， 相 对 误差 较 小 . 

2) 在 基于 一 阶 伴 随 灵敏 度 的 气动 性 能 不 确定 性 
分 析 中 ， 在 统计 意义 上 气动 偏差 统计 模型 和 外 形 偏 
差 统 计 模型 近似 满足 线性 关系 ; 对 于 满足 不 同 正 态 
分 布 的 外 形 偏差 统计 模型 气动 偏 差 概率 密度 分 布 
函数 也 近似 满足 正 态 分 布 。 

3) 基于 伴随 方法 的 由 加 工 偏差 引起 的 叶片 气动 
性 能 不 确定 性 灵敏 度 分 析 的 方法 是 有 效 的 ， 鉴 于 伴 
随 方法 在 灵敏 度 求解 方面 的 高 效 性 ， 在 进一步 提高 
伴随 灵敏 度 精 度 的 基础 上 ， 采 用 该 方法 来 预测 加 工 
叶片 气动 性 能 偏差 并 对 叶片 气动 外 形 加 工 精度 进行 
控制 优化 是 值得 期 待 的 。 
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